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Disclaimer 

Although this program was tested by its author prior to submission, 
no warranty, express or implied, is made by the author, 1620 USERS 
Group, or IBM as to the accuracy and functioning of the program and 
related program material and no responsibility is assumed by the 
author, 1620 USERS Group, or IBM in connection therewith. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter or pencil, do not use ink) 

Program No. Date 

Program Name: . 

1. Does the abstract adequately describe what the program is and what Yes No. 

it does? 

Comment ; ■ 

2. Does the program do what the abstract says ? Yes No_ 
Comment 

3. Is the Description clear, understandable, and adequate? Yes No m 

Comment \ 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

Co mment \ 

Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable? Yes No_ 

Commen t 

5. Does the source program compile satisfactorily (if applicable)? Yes No_ 

Comment^ - 

6. Does the object program run satisfactorily? Yes No_ 

Comment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No_ 

Comment • 

9. Were all necessary parts of the program received? Yes No_ 

Comment . 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt Company 
Data Corporation 

7500 Old Xenia Pike Address 
Dayton, Ohio 45432 



User Group Code, 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED 
TO PARTICIPATE IN THIS EVALUATION. 



11/09/64 



STATISTICS I 
ROBiRT J. ROBINSON 
MARQUETTE UNIVERSITY 
MILWAUKEE, WISCONSIN 
NOVSMBBR 1U, 1961 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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Program: 


STATISTICS I 




Date: 


September 21, 196l 


Programmer: 


Robert J. Robinson (Marquette University) 


Description: 


This program accepts up to £0 variables, and any number of 
observations, including missing data, and computes the standard 
deviation, standard error, coefficient of variation, number of 
observations, sum, mean, and sum of squared of deviation for 
each variable, and T-ratios and the degree of freedom between 
all variables. 

The standard deviation is computed by one of the two 
formulas, depending on the number of observations (N.) for 
the variable concerned: 1 




If N. < 30 






T-\| 






N. - 1 




if n y 30 








( *i-V 2 




The standard error is computed from: 




<5t. - - 






1 K 
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The computation of the T-ratios is optioned by setting 
a sense switch. If the computation is optioned, the T-ratio 
and the degrees of freedom are computed by one of the following 
formula sets, depending on the relation between the number of 
observations (N. and N ) of the two variables being compared: 



If N ± - H , 



10 



N.(N. -ij 
(X, - X^ \ ^ 



(£(x\ - x.) 2 ) ♦ - x.) 2 ) 

11 5 3 



D = 2(11 - 1) 



If N 
i 



T. = (X. - X.) 

1 j 1 J 



(l^N ) (K ± +-Hj -2) 



(N, + N , ) ( {£JJ - X ) 2 ) + (X.-X.) 2 ) ) 

3 3 1 J 3 



D = N 4 + N - 2 
J 



Coding 
Language: 



Fortran 



Input 
Format: 



Input is via cards. Data consists of two fixed point 
constants, N and J, which specify the number of rows in the 
input matrix (ie, the number of observations for the variables 
with the greatest number of rows), and the number of variables. 
These values are then followed by the elements of the matrix, 
punched row-wise in floating point form. Missing data is in- 
dicated by punching 9.EiiO in the proper sequence on the card. 



it 
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Thus, if the data to be used were: 



Output 
^ormat: 







y 
"2 


X 3 


\ 


1 


U.3 


6.8 


-10.9 


9.xl0 3 


2 


-6.85 


0.0 




U5.68 


3 


2.96 


U8.2 


63. 


-10.1 


H 








h.95 


5 




5.9 







The data cards could be punched as follows: 
Card 1: 



II 


j 


Xll 


X12 


X13 


XDt 


5 


h 




O 


TOT? 


9.E3 


r.2.2 


723 


J2h 


X31 


v 3 2 


X33 


0.0 


9.Eh0 


TT5758 


2.96 


1^2 


63. 


X3U 


Xhl 


Xh2 


Xk3 


Yhh 




-3 0.1 


9.2h0 


9.EU0 


9.EU0 


U.95 




X51 


7<2 


X53 


X5h 


X55 




9.E1.0 


5.9 


9.HU0 


9..ih0 


9.31)0 





721 
-6.85 



Card 2: 
Card 3: 
Card U: 



any number of elements may be punched on any card, but only the 
first 72 columns are to be used, and an element may not be split 
into two parts, (ie, cannot punch the elements Xj | ( as U.9 on the 
one card and the 5 on the next card). 

Output consists of two lines of print (or two cards) for each 
variable (ie, 2 j lines of print). The first line contains the 
following items: 

Variable number, I 
(ie, 3 if output if for 
variable X ) 



N I _ 2 



i- £_ (x T - (x T )J 

j-l i = 1 

(ie, the sum of (mean of Variable (the sum of squared 

the variable X^) t^) deviations) 



( 
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DIMENSI0NX(5O),SUM(50),SMEAN(50),SUMDV(5O) ,FNS(50) ,SUMSQ(50) 

90 READ, N, J 

N2 = N 

001 I = 1 ,J 
SUMDV( I ) = 0.0 
FNS(l)=0.0 

SUMSQ( l) = 0.0 
1SUM( I )=0.0 
3D0 12 1 = 1, J 
READ .X( I) 

IF(X(I)-9.E40)11,12,11 
1 1 FNS ( I )=FNS( I )+1 .0 
SUM(l)=SUM(l)+X(l) 

SUMSQ( I )=SUMSQ( I )+(X( I )**2) 
1 2C0NT I NUE 
N=N-1 

! F(N) 99,5,3 
5D0 6 1=1 , J 

SUMDV{ I )=(SUMSQ( I )-(SUM( I )**2/FNS( I ))) 
6 SMEAN( I )=SUM( I )/FNS( I ) 
800 10 1=1 ,J 

IF(FMS( l)-30.)17,17,19 

17 STDEV=SQR(SUMDV( I ) /( FNS( I )-1 .0) ) 
GO TO 18 

19 STDEV=SQR(SUMDV( l)/FNS( I )) 

18 STERR=STDEV/SQR(FNS( I)) 
COEFV=STDEV/SMEAN( I ) 

IF (SENSE SWITCH 2) 13,14 

13 PRINT, I ,SUM(I),SMEAN(I),SUMDV(I) 
PR I NT ,STDEV , STERR, COEFV , FNS( I ) 

GO TO 10 

14 PUNCH, I ,SUM( I ),SMEAN( I ) ,SUMDV( I ) 
PUNCH, STDEV, STERR, COEFV, FNS( I ) 

10 CONTINUE 

IF (SENSE SWITCH 1) 20,30 

20 IF (SENSE SWITCH 2) 91 ,92 

91 PRINT 
GO TO 93 

92 PUNCH 

93 J2=J-1 

DO 26 L=1 ,J2 
SMENA=SMEAN(L) 
SUMDA=SUMDV(L) 
' FNSA=FNS(L) 

21 K=L+1 

DO 25 l=K,J 

D I FF= SMENA-SMEAN( I ) 

IF(FNSA-FNS( I)) 23,22,23 

22 TRATO =DIFF*(SQR((FNS( I )*(FNS( I )-1 . 0) ) / ( SUMDA+SUMDV( I)))) 
DEGFR=2.0*(FNS( I )-1 .0) 

■ GO TO 24 

23 PART1=(FNSA*FNS( I ) )*( FNSA+FNS ( I )-2.0) 
PART2= ( FNSA+FNS( I ) )*( SUMDA+SUMDV( I ) ) 
TRATO =DIFF*(SQR(PART1/PART2)) 

b 



DEGFR=(FNSA +FNS(l)-2.0) 

24 IF (SENSE SWITCH 2) 27,28 

27 PRINT, L, I 

PRINT, TRATO, FNSA,FNS( I ), DEGFR 
GO TO 25 

28 PUNCH, L, I 

PUNCH, TRATO, FNSA,FNS( I ) , DEGFR 

25 CONTINUE 

26 CONTINUE 
30 PAUSE 

GO TO 90 
END 

END OF LISTING 
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Restrictions: 



Accuracy 
and Speed: 



Operating 
Instructions; 



The second line contains: 

or- 



(Standard Deviation 
of variable X^) 



(Standard Error 
of the mean, 
for variable Xj) 



(Coefficient 
of variation, 
for variable 



N 

I 

(number of 
significant 
items for 
variable Xj) 



If the T-test is optioned (sense switch 1 on) the following answers 
will be output when the above output is completed. 

(j-i) + (j-2) + ... + 1 sets of two lines each will be 
printed (or punched). The first line of each pair is: 



(designates 
variable Xj) 

The second line is: 
T IJ 

(T-ratio between 
and Xj) 



(desigates 

variable X ) 
J 



N I 

(number of 
significent 
items for 
variable X T ) 



(same, for 
variable 



(degrees of 
freedom between 
vari 

*.t) 



Maximum number of variables is 50. If the data is drawn from the 
same related group the formulas in this program are not valid. If the 
mean of a variable is zero, the coefficient of variation has no meaning 
as calculated. An "error E 7 11 will be printed prior to printing 
the answers for the variable whose mean is zero. All other calculations 
(including those for that variable) are valid. 

Speed not yet determined. Accuracy to 7 places. 



I. Switch Settings 

Sense Switch 1 ON Compute T-ratios. 

Sense Switch 1 OFF Do not compute T-ratios. 



Sense Switch 2 ON Output via typewriter. 
Sense Switch 2 OFF Output via card punch. 



© 
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Equipment 
Required: 



Hash 
Total: 



0,/flow switch 
Parity switch 
I/O switch 



set to program 
set to stop 
set to stop 



H. Set margins at 11 and 96. Set tabs at 28, h$, 62, 79. 
711. Put 1620 in manual mode (push reset). 

IV. Ready read hopper with the program followed by data cards. 

V. Press "Load" key. 

Program will load, and typewriter will type ' ' load data' ' . 

Merely press ' 'start' ' key on 1620 console. 

VT. When "reader no feed" light comes on, press "reader 
start" key on card reader. 



Computation will than begin, and typing will proceed. After 
the "stop" types out, if more data sets are to be run, the 
following procedure can be used: 

Press "reset" and "insert". Type h9 07500, and 
press "release" and "start". Then place data 
deck in read hopper and press "reader start". 
Proceed with the above instructions, starting with 
item VI. 

IBM 1620 card system with UOK memory. (May be used on a 
20K machine, but it will, be necessary to change the array 
sizes in the demension statement and re-compile). 

2390B827 0169927 5232$ 

(using card hash total program developed at Marquette 
Uni"ersity ). 



